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ABSTRACT 

We present a study of the effect of component separation on the recovered cosmic 
microwave background (CMB) temperature distribution, considering Gaussian and 
non-Gaussian input CMB maps. In particular, the non-Gaussian maps have been 
generated as a mixture of a Gaussian GMB map and a cosmic strings map in dif¬ 
ferent proportions. First, we extract the CMB component from simulated multifre¬ 
quency Planck data (in small patches of the sky) using the maximum-entropy method 
(MEM), Wiener filter (WF) and a method based on the subtraction of foreground 
templates plus a linear combination of frequency channels (LCFC). We then apply a 
wavelet-based method to study the Gaussianity of the recovered CMB and compare it 
with the same analysis for the input map. When the original CMB map is Gaussian 
(and assuming that point sources have been removed), we find that neither MEM nor 
WF introduce non-Gaussianity in the CMB reconstruction. Regarding the LCFC, the 
Gaussian character is also preserved provided that the appropriate combination of 
frequency channels is used. On the contrary, if the input CMB map is non-Gaussian, 
all the studied methods produce a reconstructed CMB with lower detections of non- 
Gaussianity than the original map. This effect is mainly due to the presence of in¬ 
strumental noise in the data, which clearly affects the quality of the reconstructions. 
In this case, MEM tends to produce slightly higher non-Gaussian detections in the 
reconstructed map than WF whereas the detections are lower for the LGFG. We have 
also studied the effect of point sources in the MEM reconstruction. If no attempt to 
remove point sources is performed, they clearly contaminate the GMB reconstruction, 
introducing spurious non-Gaussianity. When the brightest point sources are removed 
from the data using the Mexican Hat Wavelet, the Gaussian character of the CMB 
is preserved. However, when analysing larger regions of the sky, the variance of our 
estimators will be appreciably reduced and, in this case, we expect the point source 
residuals to introduce spurious non-Gaussianity in the CMB distribution. Therefore a 
careful subtraction (or masking) of point source emission is crucial in order to be able 
to perform Gaussian analysis of the GMB. 

Key words: methods: data analysis, techniques: image processing, cosmic microwave 
background 


1 INTRODUCTION 

One of the most valuable tools of modern cosmology is 
the observation and analysis of anisotropies in the cos¬ 
mic microwave background (CMB) radiation. A number o f 
experiments such as Boomerang iNetterfi eld et af]|2002l). 
MAXIMA jHanmij^^^J J 2 OO® , DA SI i Hrdverso^^^ll 
|20 o 3), VSA l|Grldnge je^^l200Tl . CBI ( 

ACBAR .112004. Archeops IIBenoi^^d i I 2 OO .4 

and WMAE' i^ennett et alJl2003al) have already provided 


measurements of the power spectrum of the CMB, allow¬ 
ing one to put tight constraints on the cosmological pa¬ 
rameters. In addition, the temperature distribution of the 
CMB also carries crucial information about the theory of 
structure formation. In particular, the standard inflationary 
theory predicts Gaussian fl uctuations for the CMB whereas 

S 'ogical defect models iTurok fc Snergell 1199(1 iDurrerl 
a nd non - standard inflation iLmde^^^kh^^J^U 
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I 2 OO.I iBartolo et aI]l2004^ introduce non-Gaussian signa¬ 
tures on the cosmological signal. Therefore, a detection of 
intrinsic non-Gaussianity in the GMB would have far reach¬ 
ing consequences for the current understanding of cosmol¬ 
ogy- 

When observing the microwave sky, however, one mea¬ 
sures not only the cosmological signal but a mixture of the 
GMB with other contaminant components. The most im¬ 
portant of these are the emissions from our own Galaxy 
(mainly synchrotron, free-free and dust emission), the ther¬ 
mal and kinematic Sunyaev-Zeldovich effects and the emis¬ 
sion from extragalactic point sources. In addition, the data 
will also be corrupted by instrumental noise and possibly 
some systematic effects. Therefore, in order to extract all 
the valuable information encoded in the GMB, it is critical 
to separate the cosmological signal from the other microwave 
components. This is especially relevant for the success of fu¬ 
ture CMB experiments, which will measure the microwave 
sky with unprecedented resolution, sensitivity, sky and fre¬ 
quency coverage. Most notably, these include the WMAP 
mission by NASA (that will continue to take data in the next 
few years) and the Planck mission by ESA (to be launched in 
2007), both of which produce all-sky multifrequency obser¬ 
vations of the CMB. The problem of component separation 
is a particularly important issue in the analysis of the tem¬ 
perature distribution of the CMB, since foregrounds and/or 
systematics may well introduce non-Gaussian signatures on 
the cosmological signal or, conversely, impair our ability to 
detect underlying intrinsic non-Gaussianity. The component 
separation process itself should also be well understood in 
order to avoid the introduction of artifacts that can modify 
the temperature distribution of the CMB. 

The study of the Gaussianity of the CMB has re¬ 
cently attracted a great i nterest with the re lease of 
the Ist-year WMAP data iBennett et alJ |2003^. There 


have been many analyses of the temperature distribu- 


tion 

of the CMB JKomatsu et alJ 
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and, although for some of them the data are consistent 
with Gaussian fluctuations , in other cases there has been 
a detection of non-Gaussianity and/or asymmetries, whose 
origin can not always be attributed to the presence o f fore- 
grounds or systematics. For instance, IVielva et al.l (l2004^ 
use a method based on the spherical Mexican Hat wavelet 
(MHW) and claim a non-Gaussian detection outside the 
99 per cent acceptance region, for which an intrinsic origin 
can not be discarded. By exten ding the analysi s meth od to 
orientable spherical wavelets, iMcEwen et alJ also 

obtained a significant non-Gaussian detection. 

With regard to the component separation problem, 
several methods have been proposed in the literature. There 
are basically two different approaches. The first one is to 
design a method to extract only a particular component 
from the microwave sky. The second one attempts to 
reconstruct all the components at the same time. The 


first approach is particularly well suited t o the extractio n 
of compact sources (for a review see e.g. iBarreir ^ l200fill . 
such as the emission from extragalactic point sources or 
the thermal and kinematic Sunyaev-Zeldovich (SZ) effects. 
Such m ethods include the u se of the MHW llGavon et alJ 
I 2 OOOL JVielv^^lJj20^y||H 1200,^ the matched filter 

bstal 


j|Tegm^^^^^^liveiM^os^Jl99M iHerran^^l ]l2002al. 
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2004 lMeih^B^rtie^fcM3elubroujlkj __|200^ , the scale- 
adaptive filter ISan^^Herean^^^Martfne^^TO zaled 1200 ll . 

iHerranz et al ] l20023|^^| ^^^^^^^ 1200.41 . the adap¬ 

tive top-hat filter iGhSn^^^^TiiQQ^), th e bipar ametric 
scale adaptive filter iLdpe^Ganieg^etjdj JgOOhallf ), th e 
Bayesian approach pro posed by [PtegoI^^M I 2 OO 2 II . 
the McGlean algorithm jHobson^^^^^^^^^J" 200^ or 
the wavelet-based method of IPiermol^^l ] yO^. In 
addition, there are methods that try to extract only 
the CMB componen t from the data, such a s the blind 
EM algorithm of iMartfnez-Gonzalez et ai] ll200.tl or 
the in ter nal linear rombination method jBemiettS alJ 
2003b|^_|Tegrnnxk ^e Oliveira-Costa fc Harniu^ I 2 OO 4 


Eriksen et al I l2004cll . Regarding the techniques that 


attempt to separate and reconstruct all the com- 
ponen ts at the same t i me, we have Wiener filte r 
(WF, iBoiichet et alJ Il99fil : iTegmark fc EfstathioiJ Il99fill . 
the maximum-entropy m^hod ( MEM, iHobso^^^l 
| 1998|. Il999al. IStolvarov et alJ I 2 OO 2 L 1200,4 
l20(fl ) or blind source se£axa- 


tion l|Bacci gMup^^S||200f j_|MmimetalJ |200^j_j2004 
PeUabrouiUe ^Cardos^^^a^mchoiT 1200,4 iBedinie^^ 


Patandior^^dTl200^i~ 


The effect that these techniques have on the recovered 
power spectrum of the CMB have been studied in many of 
the cases. However, no attention has been paid to the effect 
that the component separation may have on the underly¬ 
ing CMB temperature distribution. As already mentioned, 
the Gaussianity of the CMB fluctuations is predicted by 
the standard inflationary model and, therefore, a detection 
of intrinsic non-Gaussianity would have a great impact on 
our current knowledge of the universe. Thus, a complete 
understanding of the processing of the data is crucial since 
otherwise we could either misidentify spurious non-Gaussian 
signatures as true ones or impair our ability to detect intrin¬ 
sic non-Gaussianity. In the present work, we will study the 
effect of some component separation techniques on the un¬ 
derlying CMB temperature distribution. In order to do this, 
first we will apply different component separation techniques 
(MEM, WF, combined MEM-fMHW and a linear combina¬ 
tion of the frequency channels, LCFC) to simulated Planck 
data of small patches of the sky. We will perform then a 
Gaussianity analysis based on a wavelet technique to the 
reconstructed GMB map and compare the results with the 
same analysis for the input map. 

The outline of the paper is as follows. In Section we 
summarise the different component separation methods that 
are applied in this work. The wavelet technique used for the 
Gaussian analysis is explained in Sectionj^whereas Section^] 
describes the Planck simulated data. Our results are given 
in Section]^ and section [^presents some additional tests to 
understand further the role of foregrounds and noise in the 
























































































component separation process. Finally a discussion and our 
conclusions are presented in Sections |7| and |S| respectively. 


2 COMPONENT SEPARATION TECHNIQUES 

In the present section we outline the problem of component 
separation and briefly describe the techniques used for that 
purpose in this work. Our aim is to reconstruct the differ¬ 
ent components of the microwave sky, and in particular the 
CMB, in the presence of instrumental noise from multifre¬ 
quency microwave observations. 

Let us assume that we are performing a multifrequency 
observation (at n/ frequencies) of the microwave sky in a 
given direction x. We will obtain a n/-dimensional data vec¬ 
tor that contains the observed temperature fluctuations in 
this direction at each observing frequency, plus instrumen¬ 
tal noise. The measured data at the z/th frequency in the 
direction x can be written as 

ZVp Tic 

p=i 

where ric denotes the number of physical components to be 
separated and Np is the number of pixels in each map. The 
function Bi, accounts for the instrumental beam and ei/(x) 
corresponds to the instrumental noise at frequency u and 
position X. As is usual for the component separation tech¬ 
niques, we make the assumption that each of the compo¬ 
nents (except point sources) can be factorised into a spatial 
template (sp) at a reference frequency vq and a frequency 
dependence encoded in F^p. Note that point sources can 
not be factorised in this way since each of them has its own 
frequency dependence. The term rji, includes the emission 
of this component (convolved with the beam) at each fre¬ 
quency. Therefore, the component separation technique will 
attempt, at least, to extract the CMB from the data and, 
possibly, to reconstruct the rest of microwave components as 
well as a catalogue of point sources at each frequency. 

In the next subsections we outline some of the ex¬ 
isting methods to perform this component separation. We 
have considered for this work the maximum-entropy method 
(MEM), a MEM-I-MHW joint analysis, the Wiener filter 
and a method based on the linear combination of frequency 
channels (LCFC). There are other methods that have been 
pr oposed in the literature. For ins tance, the EM algorithm 
of lMartmez-Gonzalez_et_aJ. |_l|200jlMhat^_extracts only the 
CMB. and SMICA ijPellabrouiUe et alJl200,^l . which recon¬ 
struct all the components, are blind source approaches that 
first recover the power spectra of the considered compo¬ 
nents and then uses this information to form the image using 
Wiener filter. Therefore they are also, somewhat, included 
in our study. Another important technique that has been ap¬ 
plied to the pro blem of component separation in the CMB 
field is FastICA iMaino et al.ll2003l . Although this method 
is indeed very promising, some further development to im¬ 
prove the way in which FastICA deals with correlated fore¬ 
ground components, realistic beams and instrumental noise 
is, however, still needed. For this reason, we have chosen not 
to include this technique in this first analysis of the effect 
of component separation on the temperature distribution of 
the CMB. 
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2.1 The maximum-entropy method 

The maximum-entropy method (MEM) has been success¬ 
fully applied to reconstruct the microwave components 
from simulated Pla nck data in a small patch of the sky 
jHobson et alJll998l) . The algorithm was extended to deal 
with point sources (iHobson et al.l Il999a^ an d also com¬ 
bined with the Mexican hat wavelet (MHW; IVielva et alJ 
I 2 OOIIJ 1 . Subsequently, MEM has been adapt ed to work with 
spher ical data at full Planck resol ution iStolvarov et alJ 
l2002h . Finally. Istolvarov et alJ (l200f^^ have presented an im¬ 
proved MEM algorithm that can accommodate anisotropic 
noise and spatial spectral variations. All the previous al- 
gorit hms work in harmon ic (or Fourier) space. In addi¬ 
tion, iBarreiro et al.l (l2004l developed a MEM algorithm 
that combines both harmonic and real spaces and that 
can naturally deal with anisotropic noise and incomplete 
sky c overage. The method was applied to the COBE 
data. iBennett et a D (I 20031 J 1 has also applied a different 
maximum-entropy based algorithm to the WMAP data. 

In this subsection we briefly outli ne the basics of MEM . 
For a more detailed derivation, see iHobson et alJ lll99Sl . 
[T999i). If the beam is circularly symmetric and assuming 
that Fourier modes are independent, it is convenient to work 
in Fourier space (or harmonic space if dealing with all-sky 
observations). In this way, the reconstruction can be per¬ 
formed for each mode separately, what greatly simplifies the 
problem. Using matrix notation; equation 0 can be written 
for each Fourier mode as 

d = Rs -h 77 + e = Rs -h ^ (2) 

where d, rj and e are column vectors each containing n/ 
complex components and s is a column vector containing ric 
complex components. The response matrix R has dimension 
Uf xric and accounts for the effect of the beam and the spec¬ 
tral dependence of each of the components. MEM does not 
attempt to directly reconstruct the point sources but instead 
tries to prevent that this emission contaminates the rest of 
the components. Therefore, the point source contribution is 
just included in the formalism as an extra noise term. This 
is reflected in the second equality, where the instrumental 
noise and the point source emission has been combined in a 
general noise contri bution. 

As discussed in IHobson et al.l (ll99Sh . MEM should not 
itself induce correlations between the elements of the recon¬ 
structed vector. However, the microwave components may 
well be correlated. In addition, if prior information is avail¬ 
able we may wish to include it in the algorithm. In order to 
fulfill these requirements, the reconstruction is performed 
in terms of a vector of ‘hidden’ variables h. The vector of 
physical variables s is subsequently found as 


Lh , 

( 3 ) 

LL^ , 

( 4 ) 

(s(fc)s^(fe)) 

( 5 ) 


where the dagger denotes the Hermitian conjugate. C is the 
signal covariance matrix, which contains at the p-th diag¬ 
onal element, the value of the power spectrum of the p-th 
physical component at the reference frequency uo, whereas 
the off-diagonal elements contain the cross-power spectrum 
between the different components. L is a ric x nc lower trian¬ 
gular matrix, obtained by performing a Cholesky decomposi- 
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tion of C. Therefore, we can include in L our prior knowledge, 
if available, of the power spectra and cross-power spectra of 
the microwave components. Not e that L itself can b e itera¬ 
tively determined by the MEM llHobson et ahlllQQ^l . That 
is, a first reconstruction is obtained using an initial guess for 
L and subsequently the power spectra of those reconstruc¬ 
tions is estimated. These new power spectra are then used 
as a starting point for the next iteration and the process is 
repeated until convergence is achieved. 

Following Bayes’ theorem, we choose as the estimator 
h of the hidden vector to be the one that maximises the 
posterior probability given by 

Pr(h|d) oc Pr(d|h) Pr(h) (6) 

where Pr(d|h) is the likelihood of obtaining the data given 
a particular hidden vector and Pr(h) is the prior probabil¬ 
ity that contains our expectations about the hidden vector 
before acquiring any data. 

To construct the likelihood function, we assume that 
the ‘generalised’ noise contribution (containing instrumen¬ 
tal noise and emission from point sources) is well described 
by a Gaussian d istribution, which ha s been shown to work 
reasonably well iHobson et a 0119993). Therefore, the likeli¬ 
hood function is given by 

Pr(d|h) oc exp^— 

oc exp[-(d-RL/i)'''N^^(d-RL/i)] 

oc exp [-X^(h)] (7) 

where in the second line we have used 0. The noise covari¬ 
ance matrix N has n/ x n/ elements and at each fe-mode is 
given by 

N(fe) = (C(fc)C^(fc)). (8) 

Therefore, at any given Fourier mode, the vth diagonal el¬ 
ement of l\l contains the power spectrum at that mode of 
the instrumental noise plus the point source contribution to 
the z^th frequency channel. The off-diagonal terms account 
for the cross-power spectra between different channels. Note 
that if the noise is uncorrelated between channels, only the 
point sources contribute to the off-diagonal elements. 

Regarding Pr(h), MEM assumes an entropic prior prob¬ 
ability for the hidden vector h of the form 

Pr(h) oc exp[aS'(h, m)] (9) 

where ^(h, m) is the cross entropy of the complex vectors 
h and m, where m is a model vector to which h defaults in 
absence of data, a is a regularising parameter that can be 
fixed in a fully Bayesian manner. The expression for the cross 
entropy for compl ex images and the me thod for determining 
a are discussed in iHobson et al.l (Il99fjl . 

Taking into account equations 0 and 0 , maximising 
the posterior probability Pr(d|h) with respect to h is equiv¬ 
alent to minimising the function 

‘f>MBM(h) = x^(h) - aS(h, m). (10) 

Therefore the MEM reconstruction at each Fourier mode is 
given by the vector of hidden variables h that minimises the 
previous equation. 


2.2 MEM+MHW joint analysis 


As mentioned in the previous section, MEM does not di¬ 
rectly attempt to reconstruct a catalogue of point sources 
at each frequency due to the fact that this emission can not 
be factorised in a spatial template and a frequency depen¬ 
dence. Instead, point sources are included as an extra noise 
term with the aim of preventing them from contaminating 
the reconstructions of the rest of the microwave components. 
This approach has shown to work reasonably well for point 
sources with low to intermediate fluxes, greatly reducing 
the contamination due to this emission in the reconstructed 
CMB and the other microwave components iHobson et alJ 
Il999ah . However, this method does have some limitations. 
MEM does not succeed on removing the brightest point 
sources, which are present in the reconstructions although 
with significantly reduced amplitudes. This is not surpris¬ 
ing, since point source emission is modelled as an addi¬ 
tional Gaussian noise and the brightest point sources are 
not well characterised by such a model. In order to solve 
this problem, the M EM has been combined with the MHW 
iVielva et aLll^OllJl . 

The MHW technique has been developed for the extrac¬ 
tion of point source s from microwave m aps. The method was 
fir st introduced injGa yon et alJ (l200(Tll and then improved 
in IVielva et al.l ^2^^1a^ bv including the concept of opti¬ 
mal scale. The technique was then extended to deal with 
spherical data usin g the spherical Mexican hat wavelet in 
IVielva et al.l ll2nn.'jl . which provide a predicted catalogue of 
point sources from all-sky full-resolution Planck data. 

In this section we outline the basics of the MHW 
technique and how to co mbine it with the M EM (for a 
more detailed derivation see lVielva et all2001allbll . The two- 
dimensional MHW is given by 


'4j{x) = 


^/^ 1' 




(11) 


where R is the wavelet scale. By convolving the data with 
the MHW, we obtain the so-called wavelet coefficients. The 
main idea behind the MHW technique is that point sources 
are amplified in wavelet space and, therefore, they can be 
better detected. Wavelet coefficient maps can be obtained 
for each wavelet scale R. Let us assume that a point source 
of amplitude B has been convolved with a Gaussian beam of 
dispersion The value of the wavelet coefficient for scale 
R at the position of the source is given by 

{R/aaf 


w(R) 

R 


= 2y/^- 


B 


A (1 + (RK)2)2 


( 12 ) 


where A is the area occupied by the point source. 

The variance of the wavelet coefficient map at scale R 
can be calculated as 


(j^{R) = 27ri7^ 


dkP{k)k\^{Rk)\^, 


(13) 


where P{k) is the power spectrum of the map to be analysed 
and i/i corresponds to the Fourier transform of the MHW. 

As already mentioned, point sources are amplified in 
wavelet space, i.e., the level of the point sources relative to 
the dispersion of the map is higher in wavelet space than in 
real space: 


^{R) > {B/A) 
<Jni{R) CT 


( 14 ) 
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where cr is the dispersion of the observed map (dispersion in 
real space). In order to obtain the best conditions to detect 
point sources, the first term in the previous equation should 
be as high as possible. To do this, the method looks for an 
optimal scale Ropt that maximises that quantity and there¬ 
fore gives the maximum amplification of the point sources in 
wavelet space. Taking into account equations Gl and GSJ 
Ropt can be easily obtained directly from the data, since the 
optimal scale only depends on aa and the power spectrum 
of the map to be analysed. As expected, Ropt takes values 
around aa, since this is the scale that characterises the point 
sources. The background of the image also plays a role in 
determining 7?opt. Backgrounds with high power at scales 
larger than the beam size, will tend to move Ropt towards 
smaller scales than Oa and vice-versa. 

The procedure to detect point sources is as follows. 
First, the optimal scale is found and the map is filtered 
with the MHW of that scale. All pixels above a 5cru,(i?opt) 
threshold are considered as point sources candidates. The 
data map is additionally filtered with another three adja¬ 
cent scales. Those four wavelet coefficients maps are used to 
estimate the amplitude of the source by fitting the ‘experi¬ 
mental’ w(R) curve at the position of the candidate to the 
expected theoretical values given by equation I12II . If the fit 
is good, the candidates are accepted as point sources, if not, 
they are discarded. 

The combination of the MEM and MHW consists on 
applying sequentially each of the methods. First of all, the 
MHW is applied to each of the frequency maps, which al¬ 
lows one to detect the brightest point sources. These sources 
are subtracted from the original data and MEM is then ap¬ 
plied to these ‘cleaned’ maps, modelling the unsubtracted 
point sources as an additional noise term. Since the MHW 
can successfully remove the brightest sources and MEM can 
deal with those with low to intermediate flux, both meth¬ 
ods complement each other and the reconstructions of the 
different microwave components are significantly improved 
when combined. Moreover, an improved point sources cata¬ 
logue at each frequency map can also be obtained with the 
following procedure. The reconstructed maps are used to 
generate ‘mock’ data which are subtracted from the original 
ones. This provides residual maps at each frequency channel 
which contain mainly contributions from point sources and 
instrumental noise plus some residuals from other emissions. 
The MHW is then applied to these residual maps to produce 
a point source catalogue at each frequency that reaches lower 
fluxes and errors than those obtained using only the MHW. 

2.3 The Wiener filter 

Wiener filter is defined as the linear filter that min- 
imises the variance of the errors of the reconstruction (e.g. 
iR.vbicki fc Presj 11992)1 . This technique was generalised to 
deal wit h multifrequency an d mul tiresolution microwave 
data in iBoucheLeLaU tl99^ and iTegmark fc Efstathiod 
lHobsoi^^dTTi998^ derived the Wiener filter in a 
Bayesian context and compare its performance with MEM 
when reconstructing the microwave sky from simulated 
Planck observa tions, showing t hat M EM outperformed WF. 

Following iHobson et al.l (ll99Sll we will derive the 
Wiener filter within the Bayesian framework. In this con¬ 
text, the solution of the Wiener filter is found by assum¬ 


ing that the probability distribution of the sky emission is 
well described by a multivariate Gaussian characterised by a 
given covariance matrix. Therefore, the probability distribu¬ 
tion of the vector s at each Fourier mode is also described by 
a multivariate Gaussian distribution of dimension ric. This 
leads to a prior probability of the form 

Pr(h) oc exp s^Cs^ (15) 

where C is the signal covariance matrix defined in 0. Tak¬ 
ing into account the form of the likelihood found previously 
(equation GJ, the posterior probability is given by 

Pr(s|d) oc exp [—x^(s) — s^Cs] (16) 

or, equivalently, 

Pr(s|d) oc exp [—x^(h) — h^h] (17) 

where we have made used of equations (EEJ. Therefore, the 
Wiener reconstruction at each Fourier mode can be found by 
minimising, with respect to the hidden vector h, the function 

cE-wF(h) =x"(h)+hth (18) 

and then obtain the physical compone nts as s = Lh. 

We would like to point out that iHobson et al.l (Il998f) 
showed that WF can be seen as an approximation to MEM 
in the small fluctuation limit and, therefore, for Gaussian 
signals, both reconstructions become very similar. 

2.4 Linear combination of frequency channels 

If we assume that the foreground contamination is small in 
comparison with the CMB signal, a simple approach that 
combines linearly the different frequency channels (LGFC) 
in order to increase the signal-to-noise ratio can be used. In 
this case, the estimation of the cosmological signal in the 
sky at position x is given by 

Nm 

scmb(x) ^^Wj{x)dj{x) (19) 

j=i 

where Nm is the number of frequency maps to be combined 
(we do not necessarily include all the available frequency 
channels) and the data dj are given in thermodynamic tem¬ 
perature. In order to reduce the effect of instrumental noise, 
the weights of the map are chosen to be inversely propor¬ 
tional to the noise variance at each position of the sky and 
therefore the coefficients Wj(x) are given by: 


where (x) is the noise dispersion of the frequency map j at 
position X. Note that this procedure provides a map with in¬ 
creased signal-to-noise ratio and well-known noise properties 
but it does not attempt to remove any foreground contami¬ 
nation. However, when dealing with real data, some subtrac¬ 
tion of foregrounds will also be performed in the combined 
map. For instance the WMAP team has used this type of 
combination to produce a map on which performing Gaus - 
sianity studies feennett et al.[l2nn3bl : iKornatsu et alJboO.'t) 
but, before combining the maps, a fit to Galactic templates 
was removed from each of them in order to remove fore¬ 
ground contamination . In particular, they use all the re¬ 
ceivers at Q, V and W bands (at 41, 61 and 94 GHz), since 


E: 




( 20 ) 
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they are dominated, outside the Galactic plane, by the cos¬ 
mological signal. On the contrary, they avoid to include iu 
the combination the K and Ka bands (at 23 and 33 GHz) 
which contain significant Galactic emission. 

For the sake of simplicity, we will mimic the proce¬ 
dure of template fitting by reducing the Galactic foregrounds 
present in each frequency channel, prior to the linear combi¬ 
nation, down to a 10 per cent of its total amplitude. There¬ 
fore, when referring to the LCFG method, unless it is oth¬ 
erwise stated, we will assume that this foreground removal 
step has been previously performed. 


3 THE GAUSSIANITY TEST 


Given the importance of analysing the temperature dis¬ 
tribution of the CMB fluctuations, many methods have 
been developed to perform Gaussianity studies. These 
include, among others , the Minkowski functionals 

; iGott et^ |Kogu^^alJ_^9^_|Komatsue^^ 


Hea.ven.' 


tie bis» 

iSli 


iMaeueiic 


2^2^ ), th e bis pectrum Ferreira. Magneiio fc G6rskilll998t 


ST properties of hot and cold 


l2nnff) . geometrical estimators 

iBarreiro et al.l l200ll 


1200 ,1). extrema correla- 

tion function 

Naselskv & Novikov 

ll99,4lKogut et alJll99f4 

IBarreiro et al. 

Il998l iHeavens 

ShethI I 1999 II goodness of 

fit tests ijGavdi^^^ l2nn3a 


lAliapa et alJ l200,3l. 1200,11. 

multifractals JPomnilio et al. 

_ 

199,1). partition function 


analy sis lichiajag^^n~ l20nyiGhiang^N^dskv fc GolesI 

|200jJCole^e^^ 1200411 and higher criticism statistic 

In addition, wavelet techniques have been introduced 
in the last years for CMB analy ses for both small patches 
of the sky and spherical data. iFerreira. Magueiio fc ms 
ill 997^ investigate a set of stati stics based on cumulants 
and de fined in wavelet space. iHobson. .Tones fc Lasenbil 
Jl999tJl study the power of the cumulants of the dis¬ 
tribution of wavelet coefficients at each scale to detect 
cosmic strings o n small patches of the s ky. This work 
was extended by iBarreiro fc HobsonI ll2001^ who compare 
the performance of different algorithms to construct 2- 
dimensional wavelets. iForni fc Aghaniinl (ll99fl^ proposed a 
similar method that was tested o n simulated CMB maps 
containing secondary anisotropiesj[^hgnin^^Mbnii|^2^ 
and on the COBE-DMR data llAghanim. Forni fc Bouched 
I 2 OOIII . The skewness, kurtosis and scale-scale correlation co¬ 
efficients of the COBE-DMR data in the QuadCube pixeliza- 
tion using planar orthogonal wavel ets was carried out by 
iPando. Valls-Gnbaud_ fcj)an£jj[]j9^h^^ analy sis was ex- 
tended bv lMukheHee*Hobson fc Lasenbvl (1200(11 . Recently, 
wavel et methods have been ada pted to deal with sp herical 
data. IBarreiro et alJ i200flll and ICavon et al.l ll2001^ stud¬ 
ied the skewness, kurtosis and scale-scale correlation of 
the COBE-DMR data in HEALPix pixelization using the 
Spherical Haar Wavelet and the Spherical Mexican Hat 
Wavel et (SMHW), respectively. iMartfnez-Gonzalez et all 
i2002h compared the performanc e of these two sphe rical 
wavelets to detect non-Gaussianitv. lCavon et al.l ll200,3tJl put 
constraints on the /nl parameter from the COBE-DMR 
data using the SMHW. An analysis based also on the SMHW 


has b een performed on the WMAP data by IVielva et alJ 
^2004^ finding a non-Gaussian signature at scales of ~10°. 
This work has bee n extendedJjyjMuldim jjg^fc Wand ll2004ll , 
ICruz et ^ ll200!^^ and iMcEwen et ~ .l (l2005^ . Einally, steer- 
ab le wavelets on the sphere have been rece ntly proposed 
bv lWiaux. .Tacaues fc Vanderehevnstl (l2005al lbll to perform 
Gaussianity analysis of the CMB. A comparison of the per¬ 
formance of wavelet methods wit h ot her techniques can als o 
be seen in lAghanim et all 11200,311 and ICabella et aD (|200^. 

In the present work, we will use a wavelet-based 
method to de t ect no n-Gau ssianity following t h e idea s of 
iHobson et all lll999lJl and IBarreiro fc HobsorJ i200lll . In 
this section we outline the basics of the method. For a 
more detailed description of the method, see the previ¬ 
ous works. The wavelet transform has b een extensively 
described elsewhere (e.g. iDaub echiesI IT992^ . For instance, 
iBurrus. Gopinath fc Gud lll99Sll ' give a detailed introduc¬ 
tion to wavelets. The basic idea behind the wavelet trans¬ 
form is to decompose the considered signal in a series of 
wavelet coefficients that keep simultaneous information of 
real and Fourier space. Therefore each coefficient can be as¬ 
sociated to a position I and scale j of the image, what makes 
this analysis very useful in many applications. In general the 
position and scale parameters can take continuous values. 
However, for pixelized images is more couvenieut to restrict 
these parameters to take a set of discrete values. Iu partic¬ 
ular, one can construct a discrete set of wavelets that act as 
a complete basis for the digitised image. 

In 1-dimension, the wavelet basis is constructed from di¬ 
lations and translations of the mother (or analysing) wavelet 
function ip and a second related function called the father 
(or scaling) wavelet function p-. 


pj^i = 2 2^ { 2 ^ ''x — Z) 




( 21 ) 


where 2 '’ is the number of pixels of the considered discrete 
signal s{xi) and 0 ^ j ^ J — 1 and 0 ^ Z ^ 2“^ — 1 denote 
the dilation and translation indices, respectively. 

There is not a unique choice for the wavelet basis. For 
this work, we will use the so-called Daubechies-4 wavelets 
that form a re al, orthogonal au d compactly supported 
wavelet basis (see lDaubechieJl992l for the derivation of these 
functions). The discrete signal s at pixel Xi can then be writ¬ 
ten 


s{xi) = ao,opo,o{xi) + EE 


( 22 ) 


a, d correspond to the approximation and detail wavelet 
coefficients respectively. These coefficients can be obtained 
in a recursive way startiug from the data vector f{xi) = 
fflj.i, Z = 0,.., 2"^ — 1 


< 13 - 1,1 = ^ ^ h{m — 2 l)a 3 ,„ 

m 

dj-i,! = g{m - 2 l)aj,r, 


(23) 


where h, g are the low and high-pass filters associated to the 
scaling and analysing wavelet functions th rough the refine¬ 
ment equation (see e.g. lBurrus et al Idi. 

At each iteration, the vector of length 2^ is split into 
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two components: 2^“^ detail coefficients and the same num¬ 
ber of approximation coefficients. These approximation co¬ 
efficients are then used as input for the next iteration to 
construct the detail and approximation components at the 
next larger scale. The process is repeated down to the low¬ 
est resolution level considered. This leads to a number of 
wavelet coefficients equal to the original number of pixels. 
As the index j decreases from J — 1 to 0, the wavelet coef¬ 
ficients give information about the structure of the function 
on increasingly larger scales, with each scale a factor of 2 
larger than the previous one. 

The extension of the discrete wa velet transform to two- 
dimensions is not unique. Following iMallad ^1989^ we will 
construct two-dimensional bases as tensor products of the 
scaling and/or analysing wavelet at each scale j: 


142(^5 y) 

= (l>3,ii{x)4>j,iYy) 


II 

-0. 

to 


= 't’j.hOi’j.hiy) 


= i’j.hOi’rhiy) 


where H, V and D stand for horizontal, vertical and diagonal. 
The wavelet basis constructed in this way is non-redundant 
and orthogonal. 

Let us assume that we have an image with 2“^ x 2'’ pix¬ 
els. With this scheme, the image (scale J) is decomposed 
into an approximation and three detail (horizontal, vertical 
and diagonal) images corresponding to scale J —1, with each 
of the images containing 2“^“^ x 2“^“^. Note that the approx¬ 
imation image is basically a smoothed version of the input 
map, whereas the detail coefficients keep the information of 
the difference between the original and smoothed images. 
The algorithm is then applied again to the approximation 
coefficients at scale J — 1 to produce the approximation and 
detail coefficients at scale J — 2 and so on. Note that the co¬ 
efficients at scale j keep information of the structure of the 
image at scales approximately equal to the pixel size times 
~ 2-^“T 

In order to test the Gaussianity of the CMB, the 
considered temperature map is wavelet transformed using 
Daubechies-4. A certain statistic is then calculated at each 
wavelet scale for each of the three types o f wavele t coeffi¬ 
cients. In particular, following lHobson et a i] (ll999lJl . wewill 
consider the fourth order cumulant, that can be estimated 
as 


iV2[(Ar,. + l)/i4-3(iV,.-l)Ai] 

^ (Ar^._i)(iv,._2)(iV,.-3) 

A.(/,T) = {(dj-{dJ}Y). (24) 


where K 4 (/, T) corresponds to the fourth cumulant estimated 
at scale j for the T type of detail coefficient, dj are the detail 
wavelet coefficients, Nj is the number of considered coeffi¬ 
cients at scale j and T) are the estimated central mo¬ 
ments. The same procedure is repeated for a large number of 
Gaussian realisations (5000) with the same power spectrum 
as the test map. The ac 4 value obtained for the input map 
is then compared with the distribution of ac 4 obtained from 
the Gaussian realisations. Departures from this distribution 
will indicate a detection of non-Gaussianity. To avoid spuri¬ 
ous boundary effects we do not include in the analysis those 
coefficients that contain information from pixels close to the 


Frequency 

(GHz) 

Fractional 

bandwidth 

FWHM 

(arcmin) 

Pixel size 
(arcmin) 

^noise 

(xlO-®) 

30 

0.2 

33.0 

6.0 

2.0 

44 

0.2 

24.0 

6.0 

2.7 

70 

0.2 

14.0 

3.0 

4.7 

100 

0.33 

9.2 

3.0 

2.0 

143 

0.33 

7.1 

1.5 

2.2 

217 

0.33 

5.0 

1.5 

4.8 

353 

0.33 

5.0 

1.5 

14.7 

545 

0.33 

5.0 

1.5 

147 

857 

0.33 

5.0 

1.5 

6700 


Table 1. Characteristics of the Planck satellite at the 9 frequency 
channels. The central frequency and fractional bandwidth are 
given in columns 1 and 2. The antenna full-width half-maximum 
(FWHM) is given in column 3 (a Gaussian pattern is assumed 
for all channels). The characteristic pixel sizes are shown in col¬ 
umn 4. Finally, the fifth column gives the expected instrumental 
noise (assumed to be Gaussian) in AT/T per resolution element 
(which is assumed to be a square pixel of size the corresponding 
FWHM). 


borders (in particular, we do not consider ten rows/columns 
in each border). 


4 THE SIMULATED DATA 

In order to test the effect of component separation on the 
distribution of the GMB, we have generated simulated obser¬ 
vations of small patches (12.8° x 12.8°) of the microwave sky 
according to the characteristics of the Planck satellite (see 
table0. The simulations contain, in addition to the cosmo¬ 
logical signal. Galactic foregrounds (thermal dust, free-free 
and synchrotron), thermal and kinematic Sunyaev-Zeldovich 
effects, extragalactic point sources and Gaussian white noise. 

We have used two reference GMB maps for our sim¬ 
ulations. On the one hand, we have generated a Gaus¬ 
sian GMB map assuming the standard ACDM model 
(whose power spectrum w as generated using CMBFast, 
ISeliak fc Zaldarriaeal Il99(il . On the other hand, since we 
would like to test the effect of component separation on in¬ 
trinsic non-Gaussianit y, we have also used a simulated CMB 
map of cosmic strings iBouchet. Bennett fc Stebbinjll988ll . 
Given that recent CMB experiments have ruled out pure 
topological defect based scenarios fe.g. iBouchet et '^l200lll 
we will construct CMB maps with different proportions of 
the Gaussian and cosmic strings maps. Fig. Q shows the 
Gaussian CMB, the cosmic strings simulation and a map 
with a mixture of the two components. 

As Galactic foregrounds we have included simulations 
of synchrotron, thermal dust and free-free emissions, which 
have been simulated at the frequency of 353 GHz. The dif¬ 
ferent emissions have been simulated assuming that the fre¬ 
quency dependence is spatially constant. The synchrotron 
template at this frequency h as been obtained using the 
model of iGiardino et all (l2002ll . To simulate this component 
at the different Planck frequencies we have rescaled this tem¬ 
plate using a power law Iv oc with asyn = —0.9. The 

th ermal dust template has been genera ted using the model 
of iFinkbeiner. Davis fc SchlegeJ 1^9^. For simplicity, we 
have then assumed a simple grey body law with parame- 
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Figure 1. Test CMB maps used in the simulations. The different panels correspond to: Gaussian CMB (left), cosmic strings (middle), 
and Gaussian CMB plus cosmic strings in proportion 2:1 (right). The maps have been smoothed with a Gaussian beam of fwhm=5' and 
the units are thermodynamic temperature in /rK. 





Figure 2. Simulated synchrotron (left), free-free (middle) and dust (right) emissions for our reference zone (region 3) at 300 GHz. The 
maps have been smoothed with a Gaussian beam of fwhm=5' and the units are thermodynamic temperature in /rK. 


ters Td = 19>K and /3emis = 2.0 to simulate the contribution 
of this component at the Planck frequency maps. Finally, 
the free-free has been simulated using the correlation with 
dust e mission in the way proposed by iBouchet fc Gispertl 
1^9^. The template has then been extrapolated to the con¬ 
sidered frequencies using a power law with ag = —0.16. In 
order to test the robustness of the results we have tested the 
methods using four different Galactic regions of the sky with 
|fe| > 20°. The dispersion of the three Galactic components 
at the frequency of 300 GHz for the fonr considered regions 
is given in Table Two of the regions (1 and 3) are within 
the brightest 2 per cent of the sky in dust emission, whereas 
the other two regions (2 and 4) correspond to more typical 
regions with lower Galactic dust emission. Fig. shows the 
synchrotron, free-free and dust components for one of the 
considered regions of the sky. 

The kinematic and thermal SZ effects have been 
simulated using the model of IPiego et alJ ^2001^ for 
the standard AGDM model. The radio and infrared 
point sour ces have been simul ated according to the 
model of iToffolatti et alJ ll99j) (for rece n t imp rove¬ 
ments of this model, see |de_ZoHj,^ alJ I200!tI and 
Gonzalez-Nuevo. Toffolatti fc 'Auriiesn200^ for the same 
cosmological model . Finally, we have simulated the in- 



Region 1 

Region 2 

Region 3 

Region 4 

Dust 

292.9 

50.33 

115.8 

26.7 

FF 

1.46 

0.89 

0.91 

0.42 

Synch. 

1.29 

0.09 

0.28 

0.70 


Table 2. Dispersion of the three Galactic components for the 
four considered regions at 300 GHz smoothed with a Gaussian 
beam of 5 arcminutes. The units are thermodinamycal tempera¬ 
ture in fiK. 

strumental noise as Gaussian white noise according to the 
Planck characteristics. 

Planck simulated data generated from the reference 
Gaussian GMB and the Galactic foregrounds of Fig. (re¬ 
gion 3) are given in Fig.|3 Table|^gives the contribution in 
rms of each microwave component at each frequency channel 
for the same case. 


5 RESULTS 

We have performed our wavelet analysis on the recovered 
CMB maps obtained with the three considered reconstruc¬ 
tion methods. We have also studied different input CMB 
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Figure 3. Simulated Planck data containing Gaussian CMB, Galactic synchrotron, dust and free-free emissions, thermal and kinematic 
Sunyaev-Zeldovich effects, extragalactic point source emission and white noise. From left to right and top to bottom the different panels 
correspond to 30, 44, 70, 100, 143, 217, 353, 545 and 857 GHz channels. The units are thermodynamic temperature in /iK. 


Frequency (GHz) 

CMB 

kSZ 

tsz 

Dust 

Free-free 

Synchrotron 

Point sources 

30 

72.8 

0.29 

2.15 

0.25 

18.6 

29.6 

19.9 

44 

83.3 

0.37 

2.70 

0.54 

8.48 

10.3 

13.1 

70 

97.7 

0.53 

3.55 

1.45 

3.40 

2.98 

9.31 

100 

105.7 

0.67 

3.78 

3.35 

1.80 

1.23 

7.75 

143 

109.8 

0.79 

2.95 

8.35 

1.08 

0.57 

5.38 

217 

113.4 

0.95 

0.00 

31.5 

0.80 

0.31 

5.10 

353 

113.4 

0.95 

7.17 

284.8 

1.18 

0.32 

17.1 

545 

113.4 

0.95 

17.4 

5455.5 

5.10 

1.00 

287.7 

857 

113.4 

0.95 

33.2 

525918 

130.4 

18.4 

13588 


Table 3. Dispersion of the different microwave components in the Planck simulated data maps given in Fig. El (corresponding to region 
3). The units are thermodinamycal temperature in /rK. 


maps. First of all, we have produced Planck simulated data 
using our reference Gaussian CMB map and have studied if 
the reconstruction techniques preserve the Gaussian charac¬ 
ter of the temperature fluctuations (30- Secondly, we have 
considered two different non-Ganssian CMB maps (contain¬ 
ing different proportions of cosmic strings) and have studied 


if the underlying non-Gaussianity is still detected in the re¬ 
constructed temperature map To understand better 

the effect of component separation techniques on the CMB 
we have neglected the effect of point sources in these two first 
cases, since this contaminant needs to be treated in a very 
different way. In a third case, we include the emission coming 
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from point sources in the Planck simulated data and study 
the effect of this contaminant when the underlying CMB 
map is Gaussian. The study is done on reconstructed maps 
obtained using MEM and the MEM+MHW joint analysis 

In order to test the robustness of the results we have 
obtained the CMB reconstruction using simulated data of 
four different Galactic regions of the sky (see Table HJ. We 
have found that the output of the Gaussian analysis is quite 
insensitive to the choice of the Galactic region and there¬ 
fore we present our results only for simulations constructed 
using the Galactic templates of Fig. |21 We have also tested 
our results with different noise realisations, finding that the 
conclusions of the analysis are not modified. 

For the LCFC, different combinations of frequency 
maps have been tried. In particular, we have found that 
the combination of the 143 and 217 GHz channels is the 
best choice for the considered cases in the sense of retaining 
the underlying CMB temperature distribution. Therefore, 
unless it is otherwise stated, we show our results for this 
particular combination. 

5.1 Gaussian case 

The CMB reconstructed maps obtained from the Planck 
simulated data for the input Gaussian CMB are given in 
Fig. El The different panels correspond to the reconstruc¬ 
tions obtained using MEM (left), WF (middle) and the 
LCFC using the 143 and 217 GHz channels (right). These 
maps should be compared with the input Gaussian GMB 
of Fig. ID Note that the MEM and WF reconstructed maps 
look very similar. The LGFC has a lower resolution than 
the MEM and WF reconstructions, since these two methods 
perform a (partial) deconvolution of the signal. 

Fig. El shows the power spectra of the previous recon¬ 
structed CMB maps versus the input one. The solid line 
shows the power spectrum of the input Gaussian map, given 
in the left panel of Fig.0 As shown in previous works, MEM 
(dash line) can recover reasonably well the GMB power spec¬ 
trum up to ^ < 2000, whereas WF (dot-dash line) starts to 
underestimate the G^’s around I ~ 1500. The dot line cor¬ 
responds to the LCFG, where we appreciate two differences 
with respect to the input power spectrum: a defect of power 
at intermediate I's, since the combined map has a lower res¬ 
olution than the one of the input, and an excess at the higher 
multipoles, due to the presence of instrumental noise in the 
combined map. This curve can also be compared with the 
dash-three dots line, which corresponds to the power spec¬ 
trum of a combined map obtained from data channels where 
only GMB is present. This represents the best reconstruc¬ 
tion that could be achieved with the LGFC. We see that 
the reconstructed LCFC follows reasonably well this ideal 
LCFC power spectrum up to ^ ~ 1500. 

In order to test the effect of the component separa¬ 
tion technique on the underlying distribution of the CMB, 
we have applied our wavelet method to these three recon¬ 
structed CMB maps and compared the results with the in¬ 
put one. Figs. El and [7| show the value of K 4 (solid squares) 
at each wavelet scale for the input and reconstructed CMB 
maps. As already mentioned there are three different kinds 
of wavelet coefficients for the type of wavelet that we are 
using: vertical, horizontal and diagonal. Horizontal and ver- 



Figure 5. Power spectra obtained from the input and recon¬ 
structed CMB maps for the Gaussian case smoothed with a Gaus¬ 
sian beam of fwhm=5 arcmin (see text for details). 


tical details should be statistically equivalent, but we ob¬ 
tain different levels of detection because we are looking at 
a particular realisation. The i-axis on the plots gives the 
numbering k of the regions which goes as follows: region 1 
corresponds to the approximation coefficients at the lowest 
resolution; regions 2,3,4 correspond to vertical, horizontal 
and diagonal details respectively for the largest scale {j = 1); 
regions 5,6,7 give the three details in the same order for the 
next scale (j = 2) and so on. The last scale {j = 26) is not 
obtained from wavelet coefficients but corresponds to the 
value of K 4 obtained directly from the temperature map. 
Table El gives the correspondance between the region num¬ 
ber k with the scale j and type of coefficients. It also gives an 
estimation of the scale corresponding to each wavelet region. 
Since no coefficients are retained at the lowest regions, the 
x-axis runs only from fc = 8 to 26. The open circle and er¬ 
ror bars shown in the figures correspond to the average and 
the 68, 95 and 99 per cent acceptance regions obtained from 
5000 Gaussian simulations with the same power spectrum 
as the test map. For comparison purposes the distributions 
of k :4 at each scale have been normalised to unit dispersion. 
Therefore the t/-axis gives directly the values of Acr in num¬ 
bers of a (note, however, that these distributions are not 
necessarily Gaussian and thus the number of a’s in general 
will not correspond to the usual confidence intervals of the 
Gaussian case). 

Fig. El gives the wavelet statistic for the input GMB 
map (top), the MEM (middle) and the WF (bottom) recon¬ 
structions. All maps have been smoothed with a Gaussian 
beam with fwhm=5' arcmin. As expected, the value of the 
K 4 statistics for the input (Gaussian) map lie within the er¬ 
ror bars, i.e., the map is compatible with Gaussianity at all 
scales. We want to study if this behaviour is preserved in 
the reconstructed maps. We find that both MEM and WF 
produce very similar results and that the reconstructions are 
still compatible with Gaussianity. If we compare the results 
for the input and reconstructed maps, we find some small 
differences in the values of the wavelet statistics at large 
k (small scales), which are mainly due to an imperfect re¬ 
construction because of the presence of instrumental noise. 
However, this effect is small and does not modify the con- 
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Figure 4. Reconstructed CMB map obtained with MEM (left), WF (middle) and a linear combination of the 143 and 217 GHz channels 
(right) from Planck simulated data using the Gaussian GMB as input. The maps have been smoothed with a Gaussian beam of fwhm=5 
arcmin and the units are thermodynamic temperature in ^K. 


Table 4. Correspondence of region number k with wavelet scale 
j and type of coefficient. A, V, H and D stand for approxima¬ 
tion, vertical, horizontal and diagonal coefficients, respectively. 
The scale given in column 4 is calculated as 1.5'(2^“-^), where 
2-^ X 2-’ is the total number of pixels in the map. The total num¬ 
ber of coefficients and those used in the analysis for each region 
are given in the fifth and sixth columns. The rest of the coeffi¬ 
cients are discarded to avoid boundary effects. 


region 

j 

detail 

'~^scale(') 

no.coeff. 

used 

1 

1 

A 

^384 

4 

0 

2 

1 

V 

384 

4 

0 

3 

1 

H 

384 

4 

0 

4 

1 

D 

384 

4 

0 

5 

2 

V 

192 

16 

0 

6 

2 

H 

192 

16 

0 

7 

2 

D 

192 

16 

0 

8 

3 

V 

96 

64 

16 

9 

3 

H 

96 

64 

16 

10 

3 

D 

96 

64 

16 

11 

4 

V 

48 

256 

144 

12 

4 

H 

48 

256 

144 

13 

4 

D 

48 

256 

144 

14 

5 

V 

24 

1024 

784 

15 

5 

H 

24 

1024 

784 

16 

5 

D 

24 

1024 

784 

17 

6 

V 

12 

4096 

3481 

18 

6 

H 

12 

4096 

3481 

19 

6 

D 

12 

4096 

3481 

20 

7 

V 

6 

16384 

14641 

21 

7 

H 

6 

16384 

14641 

22 

7 

D 

6 

16384 

14641 

23 

8 

V 

3 

65536 

60025 

24 

8 

H 

3 

65536 

60025 

25 

8 

D 

3 

65536 

60025 


elusion that the underlying CMB distribution is consistent 
with Gaussian, according to this test. 

Fig.EI shows the results obtained from the input (top) 
and reconstructed (bottom) combined maps, constructed by 
combining the 143 and 217 GHz channels. The input com¬ 
bined map is obtained as follows: first we create two maps 
containing only CMB, filtered with the corresponding beams 
of the 143 and 217 GHz channels; we then combine the two 


maps weighting them accordingly to the noise levels of the 
two channels. This is the underlying CMB signal present in 
the reconstructed combined map and therefore we should use 
it for the comparison. Note that with this method, we can 
not aim to have better results than those obtained for the 
input combined map, since the reconstructed CMB signal 
is convolved and combined in that way. Finally, both maps 
have been smoothed with a Gaussian beam of fwhm=5 ar¬ 
cmin in order to increase the signal to noise ratio of the 
reconstructed combined map. We find also in this case that 
the input and reconstructed maps are compatible with Gaus- 
sianity. There are again some small differences between both 
plots at high k’s which are mainly due to the presence of in¬ 
strumental noise in the reconstruction. 

Therefore, if the underlying CMB map is Gaussian, all 
the considered methods preserve the Gaussian character of 
the temperature distribution. 


5.2 Non-Gaussian case 

In order to test the effect of component separation on the 
underlying CMB temperature distribution, we have also con¬ 
sidered two non-Gaussian test maps. In particular, we have 
considered mixtures of Gaussian CMB and cosmic strings 
in proportions 1:1 and 2:1 in rms. The proportion 2:1 corre¬ 
sponds approximately to the string contribution (18 per cent 
in the power spectrum) of the best fit to the CMB power 
spectru m using a mixture of inflation and cosmic strings 
found bv lBouchet et afl ||20^). 

For the non-Gaussian case with proportion 1:1, we show 
in Fig. 1^ the wavelet statistics for the 5 arcmin smoothed 
input (top), MEM (middle) and WF (bottom) maps. For the 
input map, we find very clear non-Gaussian detections for 
large k, greater than 200a at A; = 25. For the MEM recon¬ 
struction, the detections are also very clear (up to ~ 16.5 (j) 
but they have been significantly lowered with respect to the 
input map. Therefore there is a clear bias in the values of 
fi: 4 , and the CMB reconstruction tends to be more Gaussian 
than the input one. This is due to an imperfect reconstruc¬ 
tion and loss of resolution caused mainly by the presence of 
instrumental noise. At those scales where the signal-to-noise 
ratio is low, the entropic prior tends to produce a conser- 
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region number k 

Figure 6. Value of K 4 versus region number k for the case of 
Gaussian CMB. The different panels correspond to the input 
CMB (top), MEM reconstruction (middle) and WF reconstruc¬ 
tion (bottom). The open circles and the error bars correspond to 
the average and the 68 , 95 and 99 per cent acceptance regions 
obtained from 5000 Gaussian simulations with the same power 
spectrum as the test map. The solid squares are the values ob¬ 
tained for the test map. Note that the small differences in the error 
bars (of the order of about a few percent) between the input and 
reconstructions are due to the limited number of simulations. In 
addition, since the reconstruction is not perfect, the power spec¬ 
tra of the reconstructed maps differ slightly from the input one, 
what also introduces small differences in the distribution of the 
K 4 values. 



region number k 

Figure 7. Value of K 4 versus region number k for the case of 
Gaussian CMB. The panels correspond to the input (top) and 
reconstructed (bottom) combined map using the 143 and 217 GHz 
channels 


vative reconstruction with little structure, which causes the 
temperature distribution to be closer to Gaussian. A similar 
result is found for the WF reconstruction, where the highest 
non-Gaussian detection is at the level of ~ 15a. Since WF 
assumes a Gaussian prior for the signal to be recovered, this 
tends to make the reconstruction even more Gaussian than 
does the entropic prior. Fig.|^shows the results for the com¬ 
bined method. The non-Gaussianity is detected at the input 
combined map at a very high level (< 55cr). In the recon¬ 
structed combined map, the non-Gaussianity is also clearly 
detected but at a lower level (< 12a). Therefore, the re¬ 
constructed map is also more Gaussian than the input one, 
mainly due to the contamination coming from instrumental 
noise. WF and MEM produce quite similar results, although, 
in general, the MEM reconstruction produces slightly higher 
detections than WF. Also, the level of the detections ob¬ 
tained with the LCFG are lower than those of the MEM 
and WF reconstructions. Using different noise realisations 
only produces a small dispersion in the detections and does 
not modify the results. 

We have also considered map with a lower non- 
Gaussianity level constructed by mixing the Gaussian CMB 
map with the cosmic strings map in the proportion 2:1 in 
rms (right panel of Fig.0. The reconstructions obtained us- 
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Figure 10. Reconstructed CMB map obtained with MEM (left), WF (middle) and a linear combination of the 143 and 217 GHz channels 
(right) from Planck simulated data using as input a mixture of Gaussian GMB and cosmic strings in proportion 2:1. The maps have 
been smoothed with a Gaussian beam of fwhm=5 arcmin. The units are thermodynamic temperature in pK. 


ing MEM, WF and the LCFC are given in Fig.(smoothed 
with a Gaussian beam of fwhm=5 arcmin). Note that again 
the MEM and WF reconstructions are quite similar, whereas 
the LCFC has lower resolution. The corresponding recon¬ 
structed and input power spectra are also given in Fig. mi 
As before, MEM is able to recover the high £’s better than 
WF, whereas the LCFC presents an excess of power at high 
multipoles due to the presence of instrumental noise. 

Fig. II21 shows the value of /t 4 versus the region number 
for the input (top), MEM (middle) and WF (bottom) maps 
for this second non-Gaussian test map. The non-Gaussianity 
is still clearly detected in the input map although, as ex¬ 
pected, the number and level of detections are lower than 
for the mixture in proportion 1:1. In particular there are 4 
detections of non-Gaussianity outside the 99 per cent accep¬ 
tance region, with the highest detection found at fc = 25 at 
the level of ~ 50 ct. In the reconstructed maps, the level of 
the detections decrease dramatically. For MEM, we find 3 
non-Gaussian detections at level < 3.5cr, whereas for WF 
there are the same 3 detections but at slightly lower level 
(< 3 (t). Therefore there is again a large bias in the sense of 
gaussianising the reconstruction. Since we are in the limit of 
the detection, different noise realisations can make that the 
detections fall inside the 99 per cent region of acceptance. In 
any case, we find again that the MEM reconstruction pro¬ 
duces in general slightly higher detections than the one of 
WF. 

With regard to the LGFC, the results are given in 
Fig. uni Non-Gaussianity is detected at the ~ Scr level in 
the input combined map. The reconstructed map shows a 
detection outside the 99 per cent acceptance region at re¬ 
gion k = 23. However, this detection was present in the input 
only outside the 95 per cent acceptance region. This result is 
a combination of having a significant intrinsic value and 
the noise realisation and thus is not a robust detection. In 
fact, when using a different noise realisation, this detection 
disappears and other detections can appear. Therefore, we 
are again in the limit of detecting the non-Gaussianity. In 
any case, the levels of K 4 are lower than those in the MEM 
and WF reconstructions. 

Thus, in all the considered cases, there is a clear bias for 
the three studied methods and the intrinsic non-Gaussianity 
is strongly reduced in the reconstructed map. 


5.3 Effect of point sources 

For simplicity, we have not included the emission of ex- 
tragalactic point sources in the previous subsections. Now, 
we will include this contaminant in the simulations and 
study the reconstructions obtained using only MEM and 
the MEM-I-MHW joint analysis. We use the Gaussian CMB 
map given in the left panel of Fig. Q as input. 

Fig. [T31 gives the reconstructions obtained using only 
MEM (top) and the MEM-I-MHW joint analysis (bottom) 
smoothed with a Gaussian beam of fwhm=5 arcmin. The 
recovered power spectra for both cases is very similar to the 
one obtained for the MEM reconstruction of Fig. (when 
point source emission was not included). However, even if 
not seen in the power spectrum, the MEM reconstruction 
is clearly contaminated by point source emission (see top 
panel of Fig. ira . This contamination is reduced when us¬ 
ing the joint method (bottom panel of Fig. I14II . This ef¬ 
fect is also seen in Fig. EH which shows the corresponding 
wavelet statistic for the previous reconstructed maps. These 
results should be compared with the top panel of Fig. 
which shows the values of ft 4 for the input GMB map. We 
find clear detections of non-Gaussianity in the MEM recon¬ 
structed map, at the level < 60 ct. Since the underlying CMB 
signal is Gaussian, this is due to the presence of residual 
point source emission, which MEM has not been able to 
remove. However, when the MHW is used to subtract the 
brightest point sources prior to the application of MEM, this 
contamination is greatly reduced and the CMB reconstruc¬ 
tion is again compatible with Gaussianity (bottom panel of 
Fig. Similar results have been found when combining 
the MHW with WF or with the LCFC. If not subtraction 
of point sources is attempted, both methods detect residual 
point sources. In particular, for WF this detection is slightly 
lower than for MEM (since it tends to gaussianise more the 
residuals) whereas for the LCFC it is also reduced (due to 
the lower resolution of the combined map). But if the MHW 
is previously applied, the CMB reconstruction is also com¬ 
patible with Gaussianity. Therefore, the removal of point 
source emission is crucial in order to avoid the introduc¬ 
tion of spurious non-Gaussianity in the reconstructed CMB 
temperature distribution. 

As a further test, we have also applied the MEM-I-MHW 
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region number k 

Figure 8. The value of K 4 versus the region number is shown 
for the case of a mixture of Gaussian CMB and cosmic strings 
in proportion 1:1. The plots correspond to the input CMB (top), 
MEM reconstruction (middle) and WF reconstruction (bottom). 


joint analysis to Planck simulated data containing point 
sources and a non-Gaussian input CMB. In this case, we 
find that the level of the detection of non-Gaussianity is 
similar to that found in section 01 for the two considered 
proportions of cosmic strings. 
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Figure 9. The value of K 4 , for the input (top) and reconstructed 
(bottom) combined maps (143+217 GHz) for the case of a mix¬ 
ture of Gaussian CMB and cosmic strings in proportion 1:1. 



Figure 11. Power spectra obtained from the input and recon¬ 
structed CMB maps for the non-Gaussian case with a mixture 
of the Gaussian CMB and cosmic strings in proportion 2:1. The 
maps have been smoothed with a Gaussian beam of fwhm=5 ar- 
cmin. 
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region number k 

Figure 12. Values of K 4 at each wavelet region for the case of 
a mixture of Gaussian CMB and cosmic strings in proportion 
2:1: input CMB (top), MEM reconstruction (middle) and WF 
reconstruction (bottom). 


6 ROLE OF INSTRUMENTAL NOISE AND 
FOREGROUNDS 

In order to test which is the effect of each of the elements 
that take part in the component separation method, we 
present in this section some additional tests. In particular, 
we perform first some analyses to understand the effect of 
the Galactic foregrounds and the SZ emissions, without in- 



region number k 

Figure 13. The values of K 4 versus the wavelet region are 
given for the input (top) and reconstructed (bottom) combined 
(143-1-217 GHz) maps corresponding to the case of a non-Gaussian 
initial CMB generated as the mixture of Gaussian CMB and cos¬ 
mic strings in proportion 2 : 1 . 

eluding the contribution of point sources, on the CMB re¬ 
constructed distribution. Then, we consider the same case 
but reducing the level of noise by a factor of 10 and, finally, 
we further study the effect of including point sources. 

6.1 Galactic foregrounds and SZ emission 

First of all, we have applied our Gaussianity analysis to a low 
and a high frequency Planck channel (where the synchrotron 
and dust emissions are more important) without including 
point source emission in order to show how the wavelet co¬ 
efficients are affected by the presence of foregrounds. Fig. ITSl 
gives the values of K 4 for the 30 GHz (top) and 545 GHz (bot¬ 
tom) channels including the Galactic foregrounds of Fig. 
as well as our reference Gaussian CMB and the thermal and 
kinetic SZ emissions (the results are qualitatively similar for 
the other three considered Galactic regions). Since the pixel 
size at 30 GHz is 6 aremin, we consider wavelet regions only 
up to fc = 19, whereas k = 20 corresponds to the K 4 value of 
the real map. We see that, at this frequency, there are not 
deviations from Gaussianity in spite of having a certain con¬ 
tamination from synchrotron. At high Galactic latitude, this 
emission is subdominant in comparison with the CMB and 
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Figure 14. Reconstructed CMB map obtained with MEM (top) 
and MEM+MHW (bottom) for the case when point sources are 
included in the data. The reconstructions should be compared 
with the input Gaussian CMB (left panel of Fig. Q The maps 
have been smoothed with a Gaussian beam of fwhm=5 arcmin 
and the units are thermodynamic temperature in ^K. 


it is not seen in this test. We also have to point out that the 
most important contribution of the synchrotron is expected 
to occur at larger scales than those considered here, therefore 
the synchrotron could be important when analysing larger 
regions of the sky. Finally, one should bear in mind that 
the synchrotron template used in this work is based on the 
Haslam map, which has a resolution of 0.85° and therefore 
it does not contain real structure (only simulated) at the 
smallest scales, which may not be representative of the true 
emission. However, even with these uncertainties, we do not 
expect the results of this work to change with a more realistic 
simulation of synchrotron. This is so because the reconstruc¬ 
tions obtained from MEM and WF were very robust under 
variations of the Galactic foregrounds and, with regard to 
the LCFC, the synchrotron is negligible at the considered 
frequency channels (143 and 217 GHz). Regarding the 545 
GHz channel, which is completely dominated by dust emis¬ 
sion, we see very clear non-Gaussian detections at different 
scales (bottom panel of Fig. IKit . Non-Gaussian detections 
at different levels are also found at 353 GHz for the four con¬ 
sidered Galactic regions whereas at 217 GHz there is one de¬ 
tection of non-Gaussianity only for the brightest region (re¬ 
gion 1). Therefore, we see that the dust template is strongly 
non-Gaussian but for frequencies of 217 GHz or lower the 
dust emission seems to be, in general, already quite weak 
and should not affect our Gaussianity analysis (although. 
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Figure 15. Values of K 4 for each wavelet region for the 
CMB reconstruction obtained using only MEM (top) and the 
MEM-I-MHW joint analysis (bottom) when point source emis¬ 
sion is included in the data for an initial Gaussian CMB. The 
results should be compared with the top panel of Fig. which 
shows the values of K 4 for the input CMB map. 

of course, it could be more important at larger scales than 
those considered here). Of course, this is assuming that the 
uncertainties on the simulated dust emission are small. 

To further test the effect of foregrounds we have re¬ 
peated the analyses of sections 15.11 and |5 . 21 using data maps 
where the foreground components (Galactic emission and 
SZ) have been divided by a factor of 100. In this case the 
LCFC is done directly using these frequency maps (with¬ 
out further reducing the Galactic foregrounds). The results 
are very similar to those obtained in the previous section 
for both the Gaussian and non-Gaussian cases. These re¬ 
sults confirm that, for the wavelet analysis presented in this 
work and within the uncertainties of our simulations, the 
Galactic foregrounds (for high Galactic latitude regions) and 
SZ emissions do not seem to introduce any spurious non- 
Gaussianity. 

6.2 Instrumental noise 

To understand the effect of instrumental noise on the CMB 
reconstructions, we have repeated the Gaussian analyses of 
sections O and 10 reducing the noise dispersion by a fac- 
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Figure 16. Values of K .4 for each wavelet region for the 30 GHz 
(top) and the 545 GHz (bottom) frequency channels simulated 
using the Galactic foregrounds given in Fig. The 30 and 545 
frequency channels have been smoothed with a Gaussian beam of 
20 and 5 arcminutes respectively to reduce the instrumental noise. 
Note that since the 30 GHz channel has a pixel of 6 arcminutes, 
only wavelet regions with fc ^ 19 are considered, whereas k = 20 
corresponds in this plot to the value obtained from the real map. 


tor of 10 (and again without including point sources). When 
the input CMB is Gaussian, the reconstructed CMB is again 
Gaussian for MEM, WF and LCFG. The main difference is 
that the values of K 4 of the reconstructed maps are closer 
to the same values in the input map. This is especially vis¬ 
ible at the smallest scales (largest values of k) which were 
more affected by the presence of noise. With regard to the 
non-Gaussian case, the detections are clearly improved for 
the three methods and for the two proportions of cosmic 
strings considered. Fig. llTl shows the K 4 statistics for the re¬ 
constructed CMB map using MEM (top), WF (middle) and 
the combined 143-1-217 GHz map (bottom) for a mixture 
of Gaussian CMB and strings in proportion 2:1. This figure 
should be compared with Fig. ll2l and lTHl It becomes appar¬ 
ent that lowering the noise allows a much better detection of 
the intrinsic non-Gaussianity in all the cases. Therefore, the 
instrumental noise is clearly impairing our ability to detect 
intrinsic non-Gaussianity. This also indicates that MEM and 
WF perform better than the simple LCFG not only because 
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Figure 17. Values of K 4 at each wavelet region for the case of a 
mixture of Gaussian CMB and cosmic strings in proportion 2:1 
where the noise has been reduced by a factor of 10. The differ¬ 
ent panels correspond to MEM reconstruction (top), WF recon¬ 
struction (middle) and 143-1-217 GHz combined map (bottom) 
smoothed with a Gaussian beam of fwhm=5 arcminutes. 


they remove better the foregrounds but mainly because they 
attempt to denoise (and deconvolve) the signal. 
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Figure 18. Values of K 4 at each wavelet region for the 143 GHz 
channel without including (top panel) and including the emission 
of point sources (bottom panel). In order to reduce the noise, both 
maps have been smoothed with a Gaussian beam with fwhm=5 
arcminutes. 


6.3 Point sources 

In section 15.31 we have showed that point sources can be a 
very damaging foreground and that is crucial to remove, at 
least, the brightest point sources in order to avoid the pres¬ 
ence of spurious non-Gaussianity in the CMB distribution. 
As an example of the effect of extragalactic point sources, 
we show in Fig. I18l the ka statistics for the 143 GHz channel 
without including (top) and including (bottom) point source 
emission. It is clear that point sources are introducing im¬ 
portant levels of non-Gaussianity in the data, especially at 
the smallest scales, even for the 143 GHz channel where the 
contribution of point sources is expected to be relatively 
small. 

We would like to point out that we are performing 
our non-Gaussian analysis on a small patch of the sky of 
12.8° X 12.8°. However Planck will provide all-sky observa¬ 
tions of the CMB and therefore we expect to be able to 
analyse much larger regions of the sky. Since having a larger 
area of the sky will appreciably reduce the variance of the K 4 
statistics we may wonder if, in this case, we would be able 
to detect the point source residuals even after subtracting 
the brightest point sources. We have tried to give an an¬ 
swer to this question with the following qualitative argu- 
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Figure 19. Values of K 4 at each wavelet region for the 143 and 
217 GHz combined map including point source emission for a 
case where the instrumental noise has been reduced by a factor 
of 10. The combined map has been smoothed with a Gaussian 
beam with fwhm=5 arcminutes. Before combining the maps, the 
brightest point sources have been removed using the Mexican Hat 
wavelet. 


ment. For the LCFG, the part of the variance coming from 
Gaussian white noise will decrease as the considered area. 
For two-thirds of the sky this corresponds to a reduction in 
the part of the dispersion coming from the noise of the order 
of 10. Therefore, in order to mimic the case of performing 
the component separation in a larger area of the sky, we 
have constructed new data reducing the noise dispersion by 
a factor of 10 and subtracted from each channel the same 
catalogue of point sources as that of section 15.31 Then we 
have applied our non-Gaussianity test on the map obtained 
with the LCFG using the 143 and 217 GHz channel. The 
result is given in Fig. im showing a clear detection of non- 
Gaussianity. Of course, in the real case, the error bars will 
be further reduced, although not in such a straightforward 
way as for the white noise, since we are not considering the 
reduction of the variance corresponding to the GMB and the 
point sources themselves. But this means that our estima¬ 
tion can be considered as an upper limit to the error bars of 
the K 4 and thus the residual point source emission will be 
even more clearly detected. 

For the case of MEM, given that is a non-linear process, 
the way in which the dispersion of the AC 4 statistics scales 
with the area is not so simple, but we still expect our argu¬ 
ment to hold qualitatively. We have performed the Gaussian 
analysis for the data with the noise level reduced by a factor 
of 10 using MEM-kMHW and found again clear detections 
of non-Gaussianity due to the residual point sources. 

Therefore, point source emission should be carefully re¬ 
duced or masked on CMB data since otherwise it will intro¬ 
duce spurious non-Gaussianity on the analysed signal. 


7 DISCUSSION 

We have studied how different component separation tech¬ 
niques affect the underlying GMB temperature distribution. 
In particular, we have used Planck simulated data (includ¬ 
ing the cosmological signal. Galactic foregrounds, SZ effects 
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and instrumental noise) to reconstruct the CMB with MEM, 
WF and LCFC. We have also studied the effect that point 
source emission has in the reconstructed map obtained from 
MEM and a MEM+MHW joint analysis. The impact of the 
component separation method on the CMB has been quan¬ 
tified by performing a wavelet analysis of the reconstructed 
CMB map and compared it with the same analysis for the 
input map. Gaussian and non-Gaussian CMB simulations 
have been tested. We have also presented some additional 
tests to understand the effect of the different foregrounds 
and instrumental noise. 

Using as input a Gaussian CMB map (and without in¬ 
cluding the emission from point sources in the data), the 
three considered methods lead to a CMB reconstructed map 
which is compatible with Gaussianity. Therefore, none of 
them introduces spurious non-Gaussianity in the underlying 
CMB temperature distribution. However, when the input 
CMB distribution is non-Gaussian, all the methods intro¬ 
duce a strong bias in the reconstructed map, which clearly 
reduces the level of the intrinsic non-Gaussian detections 
when compared with the input map. As has been shown in 
section El this effect is mainly due to the presence of in¬ 
strumental noise which leads to an imperfect reconstruction. 

WF assumes a Gaussian prior for the signal to be re¬ 
constructed and therefore it is expected to perform well 
when the underlying GMB distribution is Gaussian. It is 
also expected that MEM will perform at least as well as WF 
since it has been shown that the entropic prio r tends to the 
Gaus sian one in the small fluctuation limit iHobson et alJ 
h Thus, it is not surprising that both methods per¬ 
form equally well when the underlying GMB is Gaussian. 
However, for non-Gaussian fields MEM would tend to per¬ 
form better than WF. This is the case for our simulations 
since we find that the level of detection of intrinsic non- 
Gaussianity is, in general, slightly higher for the MEM re¬ 
constructions than for the ones obtained with WF. Although 
for the considered cases the differences between both recon¬ 
structions are small, they could be mor e important fo r other 
type of non-Gaussianity. For insta.nce. iHobson et al ] (ll99Sll 
also show that the thermal SZ cluster profiles reconstructed 
with MEM are closer to the true ones than those obtained 
with WF. It is also interesting to note that, provided that 
we assume a good knowledge of the spectral behaviour of 
the Galactic foregrounds, their effect in the MEM and WF 
reconstructions is quite small. In fact, we have obtained our 
results for four different sets of Galactic foregrounds and 
found that the CMB reconstructions are quite insensitive to 
the input foreground maps. Regarding the LCFC technique, 
it produces lower non-Gaussian detections levels than both 
MEM and WF. We also find that, provided that we choose 
frequency channels which are relatively clean from Galactic 
foregrounds, the results do not depend on the chosen Galac¬ 
tic region (outside a Galactic cut with |6| = 20°). This result 
is confirmed by the tests performed in .section IFTTI where only 
the thermal dust seems to be playing a role at the highest 
frequency channels. However, one should bear in mind that 
we are considering only small patches of the sky and that 
the foregrounds are expected to be more important at large 
scales. In addition, our simulated templates are subject to 
uncertainties that can be especially important at subdegree 
scales. We should also point out that, although our wavelet 
technique has not been able to pick the presence of fore¬ 


ground residuals, other type of studies (e.g. phase analyses) 
may be more sensitive to the presence of foreground con¬ 
tamination. 

The fact that MEM and WF provide better reconstruc¬ 
tions than the LCFG is expected since MEM and WF are 
very powerful methods which make use of all the information 
present in the data to separate all the components. In addi¬ 
tion, they try to deconvolve and to denoise the signal, what 
clearly improves our ability to recover the CMB. However, 
they require some important assumptions which can make 
the method less robust. On the contrary, the assumptions 
required by the LCFC are less strong. When using MEM 
and WF, the frequency dependence of all the components is 
assumed to be known. This is valid for the CMB, kinetic and 
thermal SZ effects but is only approximated for the Galactic 
foregrounds. Moreover the Galactic foregrounds are assumed 
to have a pixel independent spectral index, which may be 
a reasonable approximation for small patches of the sky (as 
the ones considered here) but it does not hold when consid¬ 
ering large fractions of the sky. In addition, the results previ¬ 
ously presented have been obtained assuming full knowledge 
of the power spectra and cross-correlations between the dif¬ 
ferent components for both MEM and WF, although this 
assumption can be relaxed. 

For small patches of the sky, one could try to obtain 
directly from the data the spectral behaviour of the fore¬ 
grounds (assuming that it is spatially invariant). This could 
be done, for instance, by applying a blind-source separa¬ 
tion method. Then this information could be included in 
the MEM (or WF) algorithm to perform the final recon¬ 
struction. Another so lution is to use the ext ension of MEM 
recently developed bv IStolvarov et al.N200R|l . which success¬ 
fully accommodates spatially dust temperature variations. 
They find that the quality of the CMB reconstruction is 
comparable to the one obtained in the ideal case, when the 
dust temperature was assumed to be constant over the sky. 
The method can also be used to deal with other foregrounds. 

With regard to the influence of the assumed power spec¬ 
tra, we have repeated the analysis using MEM assuming 
neither knowledge of the power spectra nor of the cross¬ 
correlations between the components. When the input tem¬ 
perature distribution is Gaussian, the reconstructed CMB 
map is again compatible with Gaussianity. For the non- 
Gaussian cases, the level of the detections are decreased. For 
the map with Gaussianmon-Gaussian proportion of 1:1, the 
highest non-Gaussian detection goes down to ~ lOcr (versus 
~ 16.5 (t when full information was used) and for the pro¬ 
portion 2:1, there is no detection outside the 99 per cent 
acceptance region. Similar qualitative results are found for 
WF, although with lower levels of detection for the 1:1 non- 
Gaussian map (the highest detection is at the ~ 7 (j level). 
Note that these levels of detections are slightly below the 
ones obtained with the LCFC method. However, some in¬ 
formation about the power spectra of the microwave com¬ 
ponents is already available and more is expected to be ob¬ 
tained from the current and future microwave experiments. 
Therefore a realistic situation would fall somewhere between 
the two cases presented here (full information and no infor¬ 
mation). 

With regard to the LCFC, although the method is in 
principle less powerful than MEM and WF, it has the ad¬ 
vantage of not requiring so strong assumptions about the 
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underlying signals. First of all we would need to have some 
reasonably good templates for the Galactic components at 
each considered frequency channel in order to be able to 
subtract the Galactic foreground contamination down to a 
10 per cent. Given that Planck will take observations in a 
very large range of frequencies (30 to 857 GHz), including 
some channels to monitor Galactic foregrounds, this seems a 
plausible assumption. In any case, we have also repeated our 
analysis for the LCFC without subtracting any correction for 
the Galactic foregrounds and the results obtained were very 
similar to the ones presented here. Therefore, even if we were 
unable to reduce the Galactic contamination by a factor of 
10, the method would still perform similarly well. The sec¬ 
ond and most important requirement to construct the LGFC 
is a good knowledge of the properties of the noise, which is 
also a reasonable assumption for the Planck satellite. 

However, the choice of the channels to be combined is 
an important issue that must be studied in detail. For the 
present work, we have tried many different combinations of 
data maps and finally chosen to present our results with the 
combination of the 143 and 217 GHz channels. These two 
channels are relatively clean from Galactic contamination 
(which ensures that the amount of spurious non-Gaussianity 
introduced in the reconstructed CMB is very limited) and 
have a good resolution (important if we are looking for a 
non-Gaussian signal present at small scales). From all the 
tested combinations that preserved the Gaussian character 
of the CMB for the four considered regions of the sky (when 
the input was Gaussian), this was the one that yielded the 
highest non-Gaussian detections (when the input was non- 
Gaussian). We would like to remark that the optimal choice 
of channels will depend on the area of the sky that is being 
analysed as well as on the type of non-Gaussianity that we 
are trying to find. For instance, the 353 GHz channel has a 
good resolution but it may also have a high contamination 
from dust emission. However, if we are studying a particu¬ 
larly clean region of the sky, it can be helpful to add the 353 
GHz channel. Also, if we are interested in looking for non- 
Gaussian signatures at larger scales we could instead add 
other channels (e.g. the 100 GHz map). This would reduce 
the resolution of the recovered CMB map, but it will increase 
the signal-to-noise ratio at the scales of interest. In addition, 
the optimal choice may also depend on the particular test 
that is used to study the temperature distribution. 

We also have to point out that the proportion of cosmic 
strings (or other topological defec ts) present in the univ erse, 
if any, is uncertain. According to lBouchet et all (1200 ih . the 
best ht of a mixture of inflation and cosmic strings to the 
CMB power spectrum obtained from BOOMERanG (and 
other CMB experiments) yields to a string contribution of 18 
per cent (in the power spectrum). This corresponds approx¬ 
imately to a mixture of Gaussian CMB and cosmic strings 
in proportion 2 : 1 in r ms. Using the WMAP and 2dF data, 
IPogQsian et alJ ^2003^ conclude that a cosmic-string contri¬ 
bution < 10 per cent can not be ruled out. If the cosmic 
strings proportion was lower than the cases studied in this 
work, the wavelet technique would not be able to detect this 
contribution with the reconstructed methods that we have 
considered. However, note that we are using small patches of 
the sky of size 12.8° x 12.8° for the analysis. But Planck will 
provide with all-sky observations of the microwave sky and 
thus it is expected that the non-Gaussian analysis can be 


performed over a much larger fraction of the sky. Of course 
having a larger sky area implies that the proportion of cos¬ 
mic strings that can be detected is also appreciably lower. 
Therefore we expect to obtain qualitatively similar results 
to the ones found in this work for non-Gaussian maps with 
lower proportions of cosmic strings that cover larger frac¬ 
tions of the sky. 

To understand better the effect of the different contami¬ 
nants on the CMB reconstruction, the results in sections |5.1l 
and 15.21 were obtained from simulated data that did not in¬ 
clude the emission of point sources. In section 1^1^ we repeat 
the non-Gaussianity analysis including point source emission 
for an input Gaussian CMB map. In particular, we perform 
the CMB reconstruction using MEM and a MEM-fMHW 
joint analysis. We hnd that, when using only MEM, point 
source residuals contaminate the reconstructed CMB which 
clearly shows up in the wavelet statistic. However, this prob¬ 
lem is solved when using the MEM-fMHW joint analysis. 
Therefore using a tool such as the MHW to remove the 
brightest point sources prior to applying a component sep¬ 
aration technique is crucial in order to avoid that spurious 
non-Gaussianity is present in the reconstructed CMB map. 
Again, when considering a larger fraction of the sky, we may 
wonder if we would be able to see the point source residuals. 
As shown in section IQ we do expect to see the residuals 
left by this contaminant when analysing a large fraction of 
the sky, even after subtracting the brightest point sources. 
Therefore it is a critical issue to either remove or mask as 
much as possible the contribution coming from this contam¬ 
inant. 


8 CONCLUSIONS 

The aim of this work was to show the importance of consid¬ 
ering the whole component separation process when testing 
the performance of statistical tools for Gaussianity analysis. 
The method used to reconstruct the CMB map may intro¬ 
duce or erase non-Gaussianity on the cosmological signal 
and it is very important to understand this effect in order 
to give a correct interpretation of the results of the analysis. 
This type of study also allows one to establish which recon¬ 
struction technique is better, not in the sense of recovering 
the best power spectrum, but with regard to preserve the 
underlying CMB temperature distribution. 

In particular, we use a wavelet based technique to test 
the CMB recovered by three different methods: MEM, WF 
and LCFC. For the test we have used Planck simulated data 
on small patches of the sky, considering Gaussian and non- 
Gaussian CMB maps. Assuming that point sources have 
been appropriately removed, none of the methods seem to 
include spurious non-Gaussianity. However, when the un¬ 
derlying CMB is non-Gaussian, the three methods produce 
GMB reconstructions which are more Gaussian than the in¬ 
put one. This is mainly due to the presence of instrumental 
noise, which significantly affects the CMB reconstruction, 
especially at the smallest scales where the non-Gaussian sig¬ 
natures of our test maps were more important. MEM tends 
to provide slightly larger non-Gaussian detections than WF, 
whereas the detections are lower for the LCFG. This is 
mainly due to the fact that MEM and WF try to denoise 
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and deconvolve the signal and not so much to the fact that 
they remove better the foreground contamination. 

When point sources are included, if no attempt to 
remove the brightest point sources is done, they clearly 
contaminate the CMB reconstructions introducing spurious 
non-Gaussianity. However, if we hrst apply the MHW to 
remove the brightest point sources, the reconstructions are 
again compatible with Gaussianity. This shows the impor¬ 
tance of removing or masking the contamination coming 
from extragalactic point sources since otherwise it will im¬ 
pair our ability to detect intrinsic non-Gaussianity. More¬ 
over, when considering a larger region of the sky, the point 
source residuals would be more statistically signihcant and 
we expect to detect them with our wavelet technique even 
after subtracting the brightest ones. Therefore we hnd that 
this is the most damaging contaminant for Gaussian analy¬ 
sis of the CMB, at least at the scales considered in this work 
and within the uncertainties of our foreground simulations 
(that could be more important at subdegree scales). 
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